Boundary bound states in the SUSY sine-Gordon 
model with Dirichlet boundary conditions 



Chihiro Matsui 1 ^ 

1 FIRST, Aihara Innovative Mathematical Modelling Project, 
Japan Science and Technology Agency, 4-6-1, Komaba, Meguro- 
ku, Tokyo, 153-8505, Japan 

2 Institute of Industrial Science, The University of Tokyo, J^-6-l, 
Komaba, Meguro-ku, Tokyo, 153-8505, Japan 



Abstract 

We analyze the ground state structure of the supersymmetric sine- 
Gordon model via the lattice regularization. The nonlinear integral 
equations are derived for any values of the boundary parameters by the 
analytic continuation and showed three different forms depending on the 
boundary parameters. We discuss the state that each set of the nonlin- 
ear integral equations characterizes in the absence of source terms. Four 
different pictures of the ground state are found by numerically studying 
the positions of zeros in the auxiliary functions. We suggest the exis- 
tence of two classes in the SUSY sine-Gordon model, which cannot be 
mixed each other. 



1 Introduction 

Boundary systems have provided many interesting phenomena including the 
edge states and the boundary critical exponents. Also as any real materials 
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are finite systems, it is important to know how the boundaries affect on 
physical quantities. Nevertheless, analytical derivation of physical quantities 
often becomes drastically difficult due to existence of boundaries. 

Integrable systems give possibility for analytic calculation of boundary 
systems. Some examples are known to keep the system integrable even 
after imposing non-periodic boundary conditions. The famous one is the 
Heisenberg spin chain with boundary magnetic fields. This model was solved 
through the g-deformed vertex operator method [1] and the Bethe ansatz 
method [2j [3]. The interesting feature was found that the ground state 
is encoded by different structures of the roots depending on strength of 
boundary magnetic fields. 

The XXZ model with boundary magnetic fields appears also in the con- 
text of the integrable quantum field theory through the lattice regulariza- 
tion. The sine-Gordon model with the Dirichlet boundary conditions is a 
famous example, whose lattice regularized model is known as the XXZ model 
with boundary magnetic fields. As an integrable system, the S-matrix of 
the sine-Gordon model satisfies relations called the Yang-Baxter equation, 
which allows us to know the exact expressions of the scattering amplitudes. 
Also, the reflection relation is known for the integrable open system, from 
which the exact reflection amplitudes can be obtained. There are some par- 
ticles which do not scatter of each other but form bound states. For those 
cases, the procedure called the bootstrap principle has been developed to 
calculate the scattering amplitudes between bound states and their mass 
starting from the soliton mass, i.e. the mass of the lightest particle. 

However, few analytic results are available for non-periodic systems, es- 
pecially for systems with complicated symmetries such as super symmetry. 
For this reason, we study boundary effect on the ground state of the su- 
persymmetric (SUSY) sine-Gordon model with the Dirichlet boundary con- 
ditions via the lattice regularized model in this paper. Actually, knowing 
the ground state structure has the meanings more than just because the 
ground state is the most basic object. As was obtained for the spin-1/2 
case, the configuration of the ground state roots does affect on the multiple- 
integral expressions of correlation functions [H [5] . From the view point of 
the inhomogeneous XXZ model, the SUSY sine-Gordon model is the inte- 
grable extension to the spin-1 model. On the other hand, in the quantum 
field theoretic viewpoint, the SUSY sine-Gordon model is understood as the 
quantum field theory described by solitons characterized not only by the soli- 
ton charge but also by the RSOS indices. Unlike the original sine-Gordon 
case, few results are available for the SUSY case. This is partly because 
of complicated structure of the ground state, which depends on boundary 
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parameters. 

Our main aim is studying the ground state structure of the SUSY sine- 
Gordon model with the Dirichlet boundary conditions for any value of 
boundary parameters. To accomplish this aim, we first derive the nonlinear 
integral equations based on the method introduced in [6], from which one 
can read off the reflection amplitude. Some useful facts about the SUSY 
sine-Gordon model are reviewed in Section [2j In Section El the nonlinear 
integral equations are derived. Performing the analytic continuation of them 
with respect to boundary parameters, we obtain three different forms each 
of which describe scatterings on different states. After giving explanation 
about the states each set of the nonlinear integral equations is describing, 
the correct ground state is determined for full range of boundary parame- 
ters in Section [H The numerical plots of zeros in the auxiliary functions are 
given and the structure of the ground state is discussed from the analysis 
of those zeros. In the last section, we give concluding remarks and open 
problems. 



2 Supersymmetric sine-Gordon model with Dirich- 
let boundary conditions 

We first give a brief review of the SUSY sine-Gordon model. The SUSY 
sine-Gordon model is a quantum field theory described by the following 
action: 

/oo roo 
dt / dx £ S sgOM) 

1 °° ■ M M 2 (2J) 

£ssg = 2 d ^ df1 ^ + ~ ~2~ cos(/3</9)V^ + 7^2 cos2 (/^)' 

which consists of a real scalar feld ip and a Major ana fermion field tp. A 
soliton mass is denoted by M and a dimensionless coupling constant by 
f3. The value of (3 determines a main feature of the theory; the regime 
< (3 2 < Air/ 3 is known as an attractive regime and 47r/3 < /3 2 < Air a 
repulsive regime. 

Though the action (12. lj) is defined on an infinite length, it is also possible 
to consider the theory on a finite length: 

dt / dx£ SSG (x,t). (2.2) 

-oo JO 
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Special forms of boundary conditions are required in order to keep integra- 
bility of the system [7] and here we choose the Dirichlet boundary conditions: 



<p(x,t) 
<p(x,t) 



x=0 



x=L 



ip- ip(x, t) — ip(x, t) 
■-ip + ip(x,t) - $(x,t) 



x=0 



x=L 







0. 



(2.3) 



Let us remark that there is another set of the Dirichlet boundary conditions, 
although we will not consider it hereafter, given by 



ip(x,t) + i/j(x,t) 







x=0 



i/j(x,t) + ip(x,t) 



x=L 



0. 



(2.4) 



2.1 S'-matrix theory 

One of the main virtues to work on integrable systems is one can know 
the exact S-matrix by solving the characteristic relations satisfied by the 
S-matrix. The S'-matrix of the quantum field theory is defined for the 
asymptotic states: 

\AZ ai (o 1 )Aii a2 (9 2 )...Az_ ia je M )) 
= AZ ai (ei)AH a2 (e 2 ) . . . AZ"_ iau {o M )\o) 

created by applying the operators Al 3 j _ iaj (9j) to the vacuum |0). Rapidity 
of each particle 9j obeys the condition 9± > ■ • ■ > 9m for an initial state and 
9\ < • • • < 9m for a final state. The soliton creation operator has two types 
of indices one of which denotes the soliton charge o~j = ± and the other 
expresses the RSOS indices a,j = 0,^,1 under the condition \aj — aj-i\ = |. 

The S'-matrix Sssg(^) determines scattering amplitudes of a two-particle 
scattering process: 

Aim + All{0 2 ) "> Ai(9 2 ) + A^). (2.6) 

From the known properties of the S-matrix: the initial condition, the uni- 
tarity condition, the crossing symmetry, and the Yang-Baxter equation orig- 
inated from integrability of the theory together with decomposability of the 
SUSY sine-Gordon S-matrix: 



Sssg (0) = Ssg (0) ® Srsos (0), (2.7) 

the exact expressions of scattering amplitudes are obtained. We skip the 
derivation of the S-matrix, as it was closely discussed in the literatures 
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S SG (9) 



(2. 



[HJ [U Uni El H] and just show the final results. The sine-Gordon part is 
given by 

(S{6) \ 

S T (0) S r (6) 
S R {6) S t (9) 

\ 5(0)/ 

The transmissive and reflective elements of the S'-matrix St(9) and Sr(9) 
are expressed by the soliton-soliton scattering amplitude S(9) through the 
relation: 

St(9) = . m S(6>) 



Sr(9) 



sinh A(i7r — 9) 
sin7rA 



S(9) 



sinh A(i7r — 0) 
which, by using the notation u := —i9, is given by 

T(2(/ - 1)A - ^)r(2/A + 1 - v 



S(9) 



n 



Z=l 



/(« 



T((2/ - 1)A - %)T((2l - 1)A + 1 - 
On the other hand, the RSOS part is written as 

= XaMK{9) 



-u 



Srsos 
K{9) = 



a d 
b c 



1 



n 



2?ri 



)T(fc + 



2-Ki' 



The function AT^(0) is given for each set of the RSOS indices: 



X%{9) 

2 U 


= x\\(9) 

2 1 




cos 


- 


u 2 


2 


= 2 e / 2 ™cos 


(9^ 


) 


2 1 


= X%{6) 

2 u 


_ 2(«7T-S)/27ri 


cos 


f- 

\4i 


xf,{9) 

u 2 


= x\*{p) 

1 2 


= 2 e / 2 -cos 


^0 

V4i 


7T 

~ 2 



+ 



7T 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



Moreover, for the finite-size theory, scatterings on boundaries are ruled 

by the reflection matrix: 



Rssg(9)=Rsg(9)®Rrsos(9). 



(2.13) 
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The exact expression of the reflection matrix is obtained from the reflection 
relations for integrable boundaries. The sine-Gordon part of the reflection 
matrix for the Dirichlet boundary conditions is given by a 2-by-2 matrix: 



Rsg(O) - 

R ( U ) = n 



cos(£ + An) 

cos(£ — Xu] 

^)T(4A(Z - 



1=1 



r(4/A 



R (u) 

1)+1- 



cos £ 

_ 2Au\ 

7T / 



oo 



r(4/ 

cos £ 



3) A - M«r((4/ - 1)A + 1 



2Au 



cos(£ + Am) 

r(* + £ + (2Z 



1)A 



Am \ p / 1 
7T Z 1 V2 



+ (2/-l)A-^) 



r(± - | + (2/ - 2)A - ^)r(i + I + 2/A - 



/(« 



while the RSOS part is obtained as 



2 2 
OO 



P{9) = ] 



fc=i 



r(fc 



2ni 



-/(e -0) 



(2.14) 



(2.15) 



Here we gave the expression of the RSOS reflection matrix only for R1 1 (8). 

2 2 

As the ground state is non-degenerate, we assume the RSOS index for the 

ground state is |. Therefore, the other RSOS reflection factors are expected 

to be obtained for the reflection on an excited boundary. For this reason, 
i 

the expression for R^ b (0) (a,b = 0,1) will be given in the context of the 
boundary bootstrap approach in the next subsection. 



2.2 Bound states 

Particles may form bound states in the process of many-body scatterings. 
Bound states are associated with poles in the 5-matrix. The bootstrap 
approach has been developed, which allows us to know the scattering ampli- 
tudes between bound states and their masses. The same strategy is available 
to boundary reflections, which is called the boundary bootstrap. Poles of the 
reflection matrix give rapidity of boundary bound states: 



v n = - 



w N 



TT 



vr(2n + 1) 
2A 

£ vr(2iV 
A 



(2.16) 



2A 
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Only those poles which are in the physical strip Xm 9 £ (0, 5) contribute to 
boundary excitations. The condition to obtain at least one boundary bound 
state is then derived as £ > tt/2. Thus, particles trapped at the boundary 
appear when the boundary parameter exceed the threshold. Let us make a 
remark on the restriction for n and N. Through this paper, we restrict out 
interest to the repulsive regime < A < 1. Therefore, when 9 = ivq gets 
into the physical strip, the others are outside and thus we obtain only one 
boundary bound state. 

The boundary bound states form excited boundaries on which the reflec- 
tion amplitudes are read off from the boundary bootstrap. The bootstrap 
principle for the RSOS reflection factor is given by 



i) i |l) 
9\a,i\v n ) R ab( 6 ) = 9\bl\u n ) 



(2.17) 

which leads to the RSOS reflection amplitude on an excited boundary: 

Rl b (9) = P(9)K(9-iv n )K(9 + iv n )^^ [s^ cos ^ + <5 ajl _ 6 sin ^ , 

(2.18) 

where the ^-factor is the SUSY part of the boundary coupling |13|. [Hj . 



2.3 Lattice regularization 

The discretization of quantum field theories is called the lattice regularization 
and found for some models including the sine-Gordon model. Such a notion 
was first introduced by Destri and de Vega [15] . For instance, the sine- 
Gordon model can be lattice-regularized into the inhomogeneous spin-1/2 
XXZ model, or inversely, one can say that the spin-1/2 XXZ model has the 
sine-Gordon model as its low-energy effective field theory. These two models 
are equivalent in the scaling limit, where the inhomogeneity O — > 00, the 
site numbers N — > 00, and the lattice spacing a — > by keeping the soliton 
mass Ai := 4e _0 /a finite [16]. In the same way, the SUSY sine-Gordon 
model is also discretized into a lattice system. The corresponding model is 
the inhomogeneous spin-1 XXZ chain as was verified in [T7] 

Lattice regularization provided another way to analyze the continuum 
theory; though the Bethe ansatz method can be applied only to discretized 
models, the lattice sine-Gordon model admits this treatment. As relativistic 
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auxiliary space 

time direction 
A 

► space direction # ± = ± G 

Figure 1: The light-cone lattice. The inhomogeneity ±0 is assigned respec- 
tively to a right-moving particle (R) and a left-moving particle (L). 



field theory, the sine-Gordon model can be considered as a theory defined 
on the light cone, which is regarded as a Minkowski space-time. The space- 
time development on the light cone can be described by the iZ-matrices with 
inhomogeneity ±0 as rapidity of a right-moving particle (or, respectively, 
that of a left-moving particle). We skip the details which can be found in 

EE]. 

3 Nonlinear integral equations 

The scattering process on the lattice was understood through the counting 
function in the spin-1/2 case. The counting function first derived by Destri 
and de Vega |18|. [T9] gives a counting of the roots which contribute to the 
ground state, i.e. the real roots for the spin-1/2 case. With the use of 
the T-function, the object related to the transfer matrix, one can recast 
the thermodynamic Bethe ansatz equations. The sub- leading terms in the 
thermodynamic Bethe ansatz equations lead to the scattering amplitudes. 
It is known that all the elements of the scattering matrix can be obtained 
through the classification of the Bethe roots. 

Unlikely, the definition of the counting function is unclear for the spin-1 
case, due to the 2-string structure of the ground state. Instead, substitution 
has been developed in [20]. The method is based on the auxiliary functions 
[6] , which is defined through two T- functions coming from two valid transfer 
matrices T\(6) and T2(6) of the spin-1 XXZ model. 

The scattering processes can be read of from the nonlinear integral 
equations derived based on the analytic and non-zero property of the T- 
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functions. The sub- leading terms correspond to particle-scattering and 
boundary-scattering amplitudes. The source terms in the nonlinear inte- 
gral equations are understood as the effect of excitation particles. A soliton 
is obtained as an ordinary hole in the distribution of the two-strings. The 
interpretation of some other types of particles can be found in the literatures 



Such nonlinear integral equations were closely studies in [23\ YZQ [25] un- 
der the periodic boundary condition. In this section, we derive the nonlinear 
integral equations with the Dirichlet boundaries. 

3.1 Auxiliary functions 

The auxiliary functions are related to two relevant transfer matrices of the 
spin-1 XXZ model, T\{9) and T 2 (9), which are explicitly written as 



nn [22] 



T 1 (9) = h(9) + l 2 (0) 

T 2 (9) = x l (9) + x 2 (e) + x 3 (e). 



(3.1) 



Each function in T\{9) and T 2 {9) is given by 



l 2 {9) = sinh(26i - m)B_(9)(j)(9 



h{9) = sinh(26l + m)B + {9)<j){9 + m) 



in) 



Q{9 - m) 
Q{0) 

Q{9 + i7T) 

Q(0) 



(3.2) 



and 



X 1 (9) 



smh(20 - 2m)B 




) 




) 




(3.3) 
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where 



sinh M (0-e)sinh M (0 + e) 



M i 



± — ) sinh 



B±(9) = sinh 



Q(0) = \\ sinh(fl - O k ) smh(9 + O k ). 



(3.4) 



fc=i 



The variables 9k are taken as Bethe roots with positive real parts. Thus, 
the function Q{9) gives roots as its zeros. The ambiguity of 27ri-periodicity 
shall be removed by choosing the branch of a function: 



fu(x) = - In 
i 



sinh(x — iv) 
sinh(x + iv) 



as in Figure El 



(3.5) 



iv 


S > [x_ 






—iv 



TZe(x) 



Figure 2: The branch cut of the logarithm. 



We define the auxiliary functions as 
b(0) 



y(0) 



where 



Ai(0) + M9) - = X 3 (9) + X 2 (9) = 



h(0) 
T (9)T 2 (9) 

f(d) 



\i(0) 



B(O) = l + b(0) B{9) = l + b(9) Y(9) = l + y(9), 



T {9) = sinh(2#) 



(3.6) 

(3.7) 
(3.8) 

(3.9) 
(3.10) 
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It is also useful to define the ratio of h{9) and fa^d) by 

<,(*) == jf, (3.11) 

which is related to Ti(#) and 6(0) via 

T l {9) = l 2 {9){l + a^(9)) (3.12) 

6(fl) = ax (* + f ) ai (e " y) (l + ^ (* " y )) • (3-13) 

Holes in distribution of the two-strings are obtained as zeros in the transfer 
matrix T2(9) and those in T\{9) are understood as holes in real roots. On 
the other hand, the auxiliary function B(9) has zero, besides as holes, at 
two-string centers, and similarly, 1 + a\{9) at roots themselves. 

From these facts, the meanings of the auxiliary functions b(9) and a\{9) 
are understood as follows; the function b{9) gives a counting of the number 
of two-string pairs and a\{0) a counting of real roots. Indeed, the logarithms 
of the two auxiliary functions count these quantities, being associated with 
the quantum numbers Ij £ Z + \ [25] : 

i In a x (9j ) = 2nlj i In b ( 0* - — J = 2xlj . (3.14) 

As the ground state is encoded by the two-string roots, real roots should be 
considered as excited particles. For this reason, we name hole-type solutions 
of 1 + a\{9) = as type-1 holes. 

The transfer matrix T 2 (9) can be written by the auxiliary functions in 
two different ways, which we will use later to deform the contour of integrals: 

Q(9 - 2«n 

T 2 (e) = t + (9) K l\ n , I B(9) (3.15) 



Q(0 + y) 

2 

Q(0-f) 



I B(9), (3.16) 



where 



t±(0) = sinh(20±2m)S ± ffl - ^) B ± (0 + I) <A ± ^jc^ (9 ± ^ 

(3.17) 

The transfer matrix is expressed by another auxiliary function from 

the fusion relation: 

T X (9 - ^) T X (o + y)= f(0)Y(9). (3.18) 
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3.2 Classification of roots and holes 



The ground state structure of the XXZ spin chain reveals strong dependence 
on the anisotropy. For instance, when cos 7 > 1 the system is massive, while 
massless when cos 7 < 1. Here we focus on the case with t := ir/^ € 
Z>3, where the ground state is written exactly by the two-strings in the 
thermodynamic limit of the periodic system. In the context of the SUSY 
sine-Gordon model, this parameter range falls into the repulsive regime < 
A := (t — 2)" 1 < 1, in which no breathers are allowed to exist. 

Bethe roots appear either as real numbers or as complex conjugate pairs 
except for the self-conjugate roots. It is useful to classify the roots into the 
four types depending on the values of their imaginary parts: 

• Inner roots Xj (j = 1, 2, . . . , Mi) \Zm(xj)\ < 5 

• Close roots Cj (j = 1,2,..., M c ) § < \1m(cj)\ < ^ 

• Wide roots wj (j = 1,2,..., M w ) ^ < \lm( Wj )\ < 

• Self-conjugate roots w? c (j = 1, 2, . . . , Msc) \Im{w^ c )\ = 

Let us remark that pure imaginary roots or holes, which belong to neither 
of these four regimes, appear due to the Dirichlet boundaries. However, 
here we first consider the regime H± > 1 where the ground state is given 
by the pure two-string roots in spite of the Dirichlet boundaries. The exact 
two-string structure of the ground state is verified also from the numerical 
results [26j I2Z] • The reason to choose this regime is partly because we cannot 
help relying on the numerical results to know the correct ground state in 
the other regimes, in which boundary bound states may appear. Mostly, we 
show the similar statement given for the spin-1/2 case [28] that the analytic 
continuation in the Fourier transform of the nonlinear integral equations 
with respect to the boundary parameters H± naturally leads to those with 
pure imaginary objects. 

We write holes in the two valid T-matrices of the SUSY sine-Gordon 

model by h^p and hi. Since the T-matrices have zeros at the positions of 

3 J 

holes, we have the following relations: 

T 1 (hf ) ) = j = l,2,...,N 1 

' ' (3.19) 
T 2 (h j ) = j = 1,2,...,N H , 

where the number of holes and type-1 holes are given by N\ and Njj, re- 
spectively. Here we employ the assumption that both types of holes lie on 
the real axis, as in the periodic boundary system [231 EU 125] . 
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There also exist special types of holes and roots, on which B{8) crosses 
the cut of the logarithm. The imaginary part of the logarithm of b(8) has 
the derivative of a negative value: 

i\nb( Sj ) = 2nI Sj \b( Sj )\>l (i ]nb)'( Sj ) < j = l,...,N s . 

Then the logarithm of the auxiliary function B(8) should be modified as 

N S 

In 5(0) -> lnB(6) + 2m ^ H(8 - Sj ), (3.20) 

j=-N s 

where H (8) is the Heviside step function. 

3.3 Nonlinear integral equations 

Let us introduce functions T\(6) and T^iO) defined by 

U8) := T^/^Wihf}) f 2 (8) := T 2 (z)/ M (0|{fy}), (3.21) 

where fj,^(9\{h^}) has simple zeros on the real axis, besides at the origin, 

at hj and /j,(6\{hj}) at hj. Then T\{8) and T2(6) have the property of 
being analytic and nonzero near the real axis. 

Using the Cauchy theorem for the Fourier transform of T\{8) and T2(0) 
by choosing contours enclosing the real axis (Figure [3|): 

d8 e ifc ^[lnTi(0)]" = i d8 e ik ^ [In f 2 (8)}" = 0, (3.22) 



c 

the nonlinear integral equations are derived in logarithmic forms for 8 in 
the fundamental analyticity strip \Im 8\ < tt: 

In 6(0) = -(G* Cl In B)(8) - (G *c 2 In B)(8) + (d,~ fl *c 2 \nY)(6) 
+ iD hu i k (8) + iD B {8) + iD{8) + mC b 

\ny{8) = -{G [ P ] * Cl lnB)(8) + (G [ } ] * c , hiB)(0) (3 " 23) 



ITT 

^6 ~ " V 2 



+ iD SB {8) + i lim D 2 ( 8 + — - irj ) + mC y . 



Here we introduced the following notations; the convolution with integral 
contour along the real axis is written by 

/oo 
d8'f(8 - 8')g{8') (3.24) 
-oo 
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ie 


\JL 
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c 


—it 

2 





Figure 3: The integral contour of the transfer matrices T\{8) and T2(8) taken 
to enclose all the real roots and holes without any complex roots. The small 
value e is taken such that < \Xm{8j)\ — itt/2 < e, Vj. 



and for the convolution along C by 

(f*cg)(8) := / d8'f{8-8')g{8'). (3.25) 



c 

For simplicity, we express the function with a shifted variable by 

f^(8) = f(8±ir ] ). (3.26) 

The functions G{8) and G2(8) are kernels of the convolutions in the nonlinear 
integral equations: 



oc 



e 



sinh^fc l d 



G(0) := / dA; - . = 2 " = In 5(0) (3.27) 

-oo 2ir 2 smh ^-fL fc cosh | k 2m d8 

1 Z" 00 „ e" ifc ? 1 



from which one may notice that the function G{8) is related to the soliton- 
soliton scattering amplitude (|3.27p . 

The second and the fourth terms of (|3.23|) are regarded as particle source 
terms in the sense of the quantum field theory. It is useful for writing down 
the excitation particle terms D(k) and D 2 (k) to introduce effective roots: 

Oj := 8j - ysign(Xm 0,). (3.29) 
Then we have the source terms as follows: 

D(8) = z H (8) + z HW (8) - z s {8) - z c {8) - z w {8) - z sc (0) (3.30) 
D 2 (8) = (h(8) - (s(8) - Cc(0) - Cw(0) - Csc(0), (3.31) 
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where each function is given by 

A A 

z h(0)= x(0-<*j) z hW = Y X2{0-aj) 
j=-A j=-A 
A 

*s(0) = Y \X(0 - «i) + X(0 - (3.32) 

j=-A 

A A 

zc(0) = Y X(0 - a,) z w (8) = zsc(e) = Y X(9 ~ 

j=-A j=-A 

and 

A A 

CffW = £ " a i) CsW = £ ^ - aj) + X2(0 - aj)] 
j=-A j=~A 
A A 

Cc(9) = £ ^ " 5 i) CwW = &c(0) = £ X2(0 - 5,). 
j=-A j=-A 

(3.33) 

Here we denoted the number of species by A and rapidity by ay. The no- 
tation a is used for representing the complex conjugate of a. Each function 
is defined by 

X '(6) = 2ttG(9) x ' 2 (0) = 2ttG 2 (9) (3.34) 
and for the second determination of the function by 

X(0) n = X (0) + X(0 - m sign(Xm 9)). (3.35) 
The bulk term -Dbuik(^) is proportional to the number of sites N: 

smh 

D hulk (6) = N arctan — — (3.36) 
cosn H 

and each of boundary terms D-q{9) and -Dsb(^) consists of three parts: 

Db(0) = F(0; H+) + F(9; HJ) + J{9) (3.37) 
^sb(^) = F y (9; H+) + F y (9; H_) + 2G 2 U + y \ . (3.38) 

The function J{9) is given by 

f°° e - lk ^ cosh Jfcsinh^fc 

J (0) = / dfc * ^— . 3.39 

7-oo 2^ cosh^fcsinh^fc V 7 
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The functions F{6) and i^(^) have different forms depending on the values 
of boundary parameters: 



F'{9;H) 



sign(.£f) 



-Ik 



-6 



-ik 



dk 



- 



sinh - TT—^- k 



dk 



2tt 
sinh 



2 cosh Ik sinh ^k 

-w-jH+2'y ; „ 
2 ft 



2 



2vr 



2cosh^fcsinh^fc 



F'J6;H) 




i-irH N 



G 2 



2 , 



| Xm0 |<zL(idim. 



|Xm 0| < 

(3.40) 
(3.41) 



The integration constants inCb and i7rC y are determined from the asymp- 
totic behaviors of the nonlinear integral equations. 

Before ending this subsection, we give a brief comment on the nonlinear 
integral equations for 9 with the imaginary part outside the fundamental 
analyticity strip. In this case, the basic structure of the nonlinear integral 
equations is unchanged by replacing each function by its second determina- 
tion: 



In b(6) = -(G* Cl In 5X0)1! - (G *c 2 In B)(9) u + {G 2 
+ iDbuik(0)n + iD B (9) u + »£>(% + mCf 

lny(0) = -(G [ P ] * Cl ln£)(0)„ + {G [ } ] *c 2 lnS)(% 
+ i£>SB(e)n + iI>2Wn + i7rC{ 1 . 



* C2 lnF)(0)„ 



(3.42) 



The integral constants iTrCjf and mC^ are determined in a similar way to 
i7rCfe and mC y from the asymptotic behaviors, though in general, they have 
different values for these two cases. 



3.4 Sum rules 

We discussed so far the nonlinear integral equations without referring the 
restriction of the number of roots and holes. Actually, they are not arbitrary 
but are constrained by the so-called sum rules. The sum rules are derived 
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from the analysis of asymptotic behaviors of the nonlinear integral equations: 



N H - 2N S = M c + 2{S + M W + M SC ) 



3 7 (2S + # + l) 7(2S + F + 1 



TT 



TT 



Ni - 2N§ = S + M W + M sc + - M R 



j(2S + H + l) 



TT 



(3.43) 
(3.44) 



where N$ is the number of real special objects and M c is the number 
of close roots whose imaginary parts lie on (tt/2,tt) or (— tt, — tt/2). The 
notation |_^J was used to express the largest integer part of x. 

3.5 Scaling limit 

The nonlinear integral equations of the lattice SUSY sine-Gordon model are 
identical to those for the original SUSY sine-Gordon model in the scaling 
limit. The scaling limit is realized by taking the three limits N, O — > oo and 
a — > (aN := L) by keeping the soliton mass Ai finite. These quantities 
appear only through the bulk term (13.23P which results in 



./V arctan 



sinh# 

COsh scaling limit 



2iMLsinh9. 



(3.45) 



4 Boundary bound states and the ground state of 
BSSG model 

In the context of the S'-matrix theory, rapidity of particles bounded at the 
boundaries emerge as poles in the reflection matrix. On the other hand, 
boundary bound states are obtained as pure imaginary roots from the Bethe 
ansatz context. These two objects nevertheless does not coincides in the 
finite system, due to the finite-size deviation of the Bethe roots. They have 
the same values only in the infinite- volume limit ML — > oo, which we are 
going to consider in this section. 



4.1 Large volume limit 

In the relation with the quantum field theory, we call the ML — > oo limit the 
infrared (IR) limit. The nonlinear integral equations (13. 23ft are simplified in 
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the IR limit into 



In b(6) = (Ga 2 *c 2 \nY)(9) + 2iMLsmhO + iD B (0) +iD(0) +mC b 
hxy(9) = iD SB (9) + iD 2 {9) + mC y . 

(4.1) 

One may notice that we also took the scaling limit. 

Suppose the Dirac sea consists of the pure two-strings without particles 
for H± > 1. The one-soliton excitation is then described by appearance 
of a hole without any other objects; Then the nonlinear integral equations 
with a hole are equations for the behavior of one soliton in a finite interval 
L. Therefore the sub-leading terms of the nonlinear integral equations (|4.ip 
give the scattering amplitude on the ground state boundary: 



exp[(G^ * l * C2 lnY)(6)+tD B (6)} 



H->1 S{29) 

(4.2) 

from which one can read off the relations between boundary parameters of 
the S-matrix theory and the lattice-regularized model: 

t-H = l + ^- 

1 nX (4.3) 

Here we set a superscript \B) <g> \a) to 7£ssg> whose first component indicates 
the sine-Gordon part of a boundary state and the second the RSOS part. 
There are no solitons as the state |0) represents with the RSOS index ^ due 
to non-degeneracy of the ground state [29]. We will denote the existence of 
a pure imaginary root 9 = ivq by a superscript and the corresponding 
RSOS index by one of the doublet RSOS states \a) (a = 0, 1). 

From now on, we call the nonlinear integral equations without the source 
terms as the bare nonlinear integral equations. Once again the bare non- 
linear integral equations given in (14.ip were written for the pure two-string 
ground state. On the other hand, the parameter relations (|4.3j) indicate that 
the first pole in the reflection matrix 6 = wq gets into the physical strip for 
< H- < 1, that is, a boundary bound state emerges in this regime. 

Now the question is which state the bare nonlinear integral equations 
are describing. For simplicity, we let only H_ move by fixing H + at a 
value larger than 1. Thus, boundary bound states would appear only at 
the (— )-boundary. Let us first sum up three different forms of the boundary 
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(a) 


H > 1 


Db(0)\h>i 


Dsb(0)\h>i 


(b) 


1 >H > -1 


Db(6)\i>h>-i 


D S b(0)\i>h 


(c) 


-l> H 


D b (0)\-x >h 


Dsb{8)\i>h 



Table 1: Three regimes with different forms of boundary terms. 



term Db{0) obtained by analytic continuation with respect to iT_ (Table [T]). 
In the case of (a), we obtained (|4.ip describes the pure two-string ground 
state. In the regime (b), the analytic continuation with respect to H- makes 
Db{0) into a different form. Recall the boundary bootstrap equation for the 
sine- Gordon part: 

K^(0;\,£) = K^(9;\,£)S(9-iv )S(9 + iv ), (4.4) 

which precisely gives the expression of Db{9) for the regime (b). This means 
the analytic continuation of the nonlinear integral equations naturally in- 
volves the appearance of boundary excitation as the appearance of a hole at 
9 = wq. 

Thus the straightforward computation leads us to obtain that the re- 
flection amplitude on the first excited boundary, given by the appearance of 
pure imaginary holes at 9 = ±7r(l — H^)/2, coincides with what is described 
by the bare nonlinear integral equations for 1 > > — 1: 



exp[(G^ 2l *c 2 ^Y)(9)+iD B (9)} 



n^ a \9 ] x,^)n^ ) (9-,x,^) 



i>ff_>-i S(29) 

(4.5) 

On the other hand, the bare nonlinear integral equations for i?_ < — 1 
are related to those for H_ > 1 by 



(Gp l * C2 \nY)(9) + D B (9) 

— 1>X1 

(dr fl *c 2 lnY)(9) + D B (9) 



+ X2(9 — )+X2[9 + 



H->1 V 2 / V 2 

(4.6) 

which implies the appearance of type-1 holes at 9 = —iirH_/2. As was 
already referred, a SUSY sine-Gordon particle is interpreted as a hole in the 
2-strings. Since the type-1 hole lies in the distribution of real roots, this 
should be considered as an excitation for the system with the ground state 
described by the pure two-strings. 
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We give an interpretation on this situation; the pure two-string ground 
state or states consisting of ordinary holes are realized by an even number 
of Bethe roots. If the system is encoded by an odd number of Bethe roots 
from the beginning, not every root can make pairs and thus strong enough 
boundary field arrests only one root. Thus, our suggestion is that the bare 
nonlinear integral equations for iJ_ > 1 and i?_ < — 1 describe difference 
sectors of the SUSY sine-Gordon model, which cannot be mixed simply by 
parameter transformations or particle excitations Q 

4.2 Ground state 

In this section, we discuss the ground state of the SUSY sine-Gordon model 
with Dirichlet boundary conditions on a finite interval. The similar analysis 
has done for the spin- 1/2 chain under existence of boundary magnetic fields 
|30| and they found 10 regimes with different structures for the ground state. 

The main protagonists of the analysis are pure imaginary roots/holes 
and poles in the auxiliary functions. As is referred in the previous section, 
the poles which may appear for a given value of H± are naturally taken into 
consideration by analytic continuation of the nonlinear integral equations. 

The degree analysis indicates that being set as w = e e , there exist N +1 
holes in 1 + a\{9) and N + 2 holes in B{9) for the system of length N under 
existence of N roots, as is also numerically checked in Figure HI Now we look 
into how these holes move depending on the values of boundary parameters. 
Following the previous section, we let only one of the boundary parameters 
H- move by fixing the other H + at arbitrary value larger than 1. For 
boundary field H- larger than 1, all holes lie on the real axis (Figure [4.2p . 
As boundary field strengthens, that is, as the value of H- becomes smaller, 
one of the holes in T2(6) moves toward the origin and then on the imaginary 
axis after H- passes through 1 from above (Figure 14. 2p . This hole seems 
to be fixed at 6 = m{l — H_)/2 up to exponentially small correction of the 
system size. At the same time, one of the poles in the auxiliary function 
B{6) locating at 6 = — wr(l — JT_)/2 also passes the origin and subsequently 
the Bethe equations admit a pure imaginary root. The appearance of a pure 
imaginary root is understood in the context of the Bethe ansatz equations 
for the large enough system as the pole with a positive value of the imaginary 
part is compensated by the term of the TVth power which vanishes as N gets 
large. As a result, the Bethe ansatz equations with large enough N admit 

"The bare nonlinear integral equations for H- > 1 cannot reproduce all the eigenvalues 
of the RSOS part of the scattering matrix [24| . We expect the two of them which did not 
show up would be obtained from those for H < — 1. 
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pure imaginary roots for H_ < 1. Here we remark that a hole lies on the 
distribution of the two-strings. Thus, it should considered that two roots 
are lost when a hole emerges. Assuming the value of hole rapidity is given 
by string-ceter rapidity of two roots, we conclude that two imaginary roots 
appear at 8 = —iirH_/2 and —m(H_ — 2)/2. One can easily check that the 
wave function composed by imaginary rapidity exponentially small from the 
boundary. For this reason, the imaginary roots can be understood to form 
boundary bound states. 

Now we consider smaller H- than 1. When crosses 0, one of the 
holes in Ti(9) reaches the origin moving on the real axis. This hole then 
moves up along the imaginary axis as H- decreases fixed at 9 = —iTrH_/2 
up to exponentially small finite-size correction (Figure \4.2\i . The auxiliary 
function 1 + ai(8) has a pole of the same value as the hole-type zero of 
T\{9) up to the correction and picks it up to the imaginary axis. Thus 
a pure imaginary root 9 = —inH/2 appears. One has to be careful that 
this time the hole emerges from 1 + a\(9), which counts real roots. Thus, 
the appearance of a type-1 hole is achieved simply by removing one root. 
Therefore, a type-1 hole generates a different boundary bound state from 
the two-imaginary-root case. 

Finally when H_ becomes smaller than —1, the other hole of T2{9) 
reaches the origin and goes up on the imaginary axis (Figure 14. 2p . Being 
fixed at 9 = — wr(-ff_ + 1), this hole is picked up by a pole in the auxiliary 
function B(9). Thus another imaginary roots appear at 9 = —iirH_/2 and 
—iir(H_ + 2)/2 in this regime. 

4.3 Symmetry with respect to boundary parameters 

We give some comments on the symmetries of the SUSY sine-Gordon model. 
From the action of the model ()2.2|) . we obtain that the theory is invariant 
under the transformations: 



In comparison with the isotropic limit of the model, it is useful to rewrite 
the symmetry in terms of the lattice boundary parameter H± by utilizing 
the relations (14. 3D : 



(a) <p± -> 27r//3-^ ± e± -> 2£o-£± H± -> -H±-2 



2-ir 
2vr 



(4.7) 



(4.8) 




<P± ->• ~¥>±- 



(4.9) 
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me) 



-ID -05 



(a) (left) Zeros of l+ai(0); (right) B{6) for H+ = 1.5 and H- = 2.2. 



me) 



(b) (left) Zeros of l+ai(0); (right) B(0) for #+ = 1.5 and if_ = 0.3. 



mi) 



mi) 



■ ■ -05 - • - 

-10 - 

(c) (left) Zeros of 1 + ai(0); (right) B(6») for = 1.5 and H- = 
-0.5. 



. m) . . . u 

ID 


- . . . . . me) . . . , 




-05 


■ -05 
-ID 


05 


-05 


05 



(d) (left) Zeros of 1 + ai(0); (right) B[9) for H+ = 1.5 and H- = 
-1.8. 



Figure 4: Zeros of the auxiliary functions 1 + a\{9) and B(9) is plotted 
for N = 8 with 8 roots in the homogeneous and isotropic limit. Roots are 
plotted by green dots and hole-type zeros by red dots. 
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(b) <^± 

(c) V± 



2tt//3 + ip ± 



-ip± 



e± -»• 2£o + e± H± ^ H ± - 2t 

e± -> -£± #± -> -#± + 2(t - 1) 



where 2£ := vr(2A + 1) and £± = 2mp±//3. 

The corresponding spin system to the SUSY sine-Gordon model has the 
following interactions; 



H 



N-l 
5=1 



.7 J 



-i + Wb- 



The bulk term is given by 



=T j - (Tj) - 2(sin 7 )^ [T? + (S*y + (S 



'j+V 



z\2] 



7 



where 



rp Q Q rpj_ _ qX qX . qy qy 

3 — 3 ' J+l 3 — j J+l ' j 7+1 



rpZ 



qZ qZ 



and the boundary term by 

n B = h 1 (H^S z 1 +h 2 (H^)(S z 1 ) 2 



where 



h 2 (H) = I sm(2 7 ) 



cot 



2 

cot ( 



h 1 (H + )S' N + h 2 (H+)(S z N ) 2 , 



V 2 



(f) + co t , 



V 2 



+ 7 
+ 7 



(4.10) 



(4.11) 

(4.12) 
(4.13) 

(4.14) 
(4.15) 



In Figure El the boundary energies of the spin chain are depicted. One ob- 
tains the same symmetries of the SUSY sine-Gordon model in the boundary 
energies of the spin chain. Thus, the three symmetries of the quantum field 
theory survive even after taking the isotropic limit, which supports validity 
of our numerical results in the analysis of the SUSY sine-Gordon model. 



5 Concluding remarks 

In this paper, we studied the ground state of the SUSY sine-Gordon model. 
The structure of the ground state roots strongly depends on the values of 
boundary parameters. What was found here is that there are four regimes 
which are characterized by different features of Bethe roots. The detailed 
structure is; 
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Figure 5: Energy of boundary terms in the spin-1 XXZ model with boundary 
magnetic fields. 



2. in the regime either H- > 1 or < —2, the ground state is described 
by the pure two-strings. 

3. For 1 > H— > 0, the ground state include one soliton with rapidity 
9 = 27r(l — HJ)/2. In the context of the Bethe ansatz equations, the 
ground state is encoded by the two-strings and two imaginary roots 
9 = -inH^/2 and -m{H_ - 2)/2. 

4. In the regime > H— > — 1, a type-1 hole contributes to the ground 
state. Rapidity of this hole is given by the corresponding Bethe root 
itself 9 = —mH^/2. This situation is realized only when the original 
system is characterized by the odd number of roots. Thus, the lowest 
energy state in this regime belongs to a different class from the other 
regimes. 

5. Finally for — 1 > H— > —2, the ground state is given again by the 
two-strings and one soliton with pure imaginary rapidity, but at 9 = 
—m{H- + l)/2 this time. Correspondingly, the ground state includes 
two pure imaginary roots at 9 = —iirH_/2 and —m(H- + 2)/2. 

We also discussed what is described by the analytic continuation of the 
nonlinear integral equations with respect to the boundary parameters in the 
Fourier space. It was found that 

1. the bare nonlinear integral equations for H > 1 describes the two- 
string ground state. 
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2. In the regime 1 > H_ > —1, they are written for the state which 
includes a hole at 9 = m{l — H_)/2. 

3. the bare nonlinear integral equations for < — 1 are the equations 
for the state with a type-1 hole at = —iirH^/2, which belongs to the 
different class of the above two. 

The nonlinear integral equations provide the sum rules, the restriction 
on the number of various roots and holes, from their asymptotic behavior. 
The sum rules help us to know the possible configurations of roots and holes. 
We expect that the all particle and reflection amplitudes of the SUSY sine- 
Gordon model would be recasted from configurations of holes and roots, 
which include close roots or wide roots. It is also an interesting open problem 
to study scatterings on the system whose ground state is encoded by the odd 
number of roots, which may give the lacked two eigenvalues of the RSOS 
part of the reflection matrix. 

We briefly remarked the symmetry of the SUSY sine-Gordon model with 
respect to the boundary parameters. They are translated into the symme- 
tries of lattice boundary parameters via the identification obtained from the 
comparison of the sub-leading term of the nonlinear integral equations with 
the reflection amplitudes. As it was found that these symmetries survive in 
the homogeneous limit from the energy plots of the boundary terms, we ex- 
pect the ground state of the spin chain has the same structure as the SUSY 
sine-Gordon model. 
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